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Abstract 

We propose a generalized feedback model and compressive sensing based opportunistic feedback 
schemes for feedback resource reduction in MIMO Broadcast Channels under the assumption that both 
uplink and downlink channels undergo block Rayleigh fading. Feedback resources are shared and are 
opportunistically accessed by users who are strong, i.e. users whose channel quality information is above a 
certain fixed threshold. Strong users send the same feedback information on all shared channels. They are 
identified by the base station via compressive sensing. Both analog and digital feedbacks are considered. 
The proposed analog & digital opportunistic feedback schemes are shown to achieve the same sum-rate 
throughput as that achieved by dedicated feedback schemes, but with feedback channels growing only 
logarithmically with number of users. Moreover, there is also a reduction in the feedback load. In the 
analog feedback case, we show that the proposed scheme reduces the feedback noise which eventually 
results in better throughput, whereas in the digital feedback case the proposed scheme in a noisy scenario 
achieves almost the throughput obtained in a noiseless dedicated feedback scenario. We also show that 
for a given fixed budget of feedback bits, there exists a trade-off between the number of shared channels 
and thresholds accuracy of the fed back SNR. 

Index Terms 

Compressed sensing, feedback, lasso, multiple-input multiple-output (MIMO) systems, opportunistic, 
protocols, scheduling. 

I. Introduction 

Recently, it has been shown that dirty paper coding (DPC) achieves the sum-rate throughput of the 
multiple-input multiple-output (MIMO) broadcast channel fl], El- However, it requires a great deal of 
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feedback as the transmitter needs perfect channel state information for all users and is computationally 
expensive 0. Since then, many works have attempted to achieve the same sum-rate throughput with 
imperfect channel state information (reduced feedback load). This was done by applying opportunistic 
communication in the forward link Q- ifTOil . By reviewing the feedback protocols suggested in literature, 
one can note that generally the following three components are fed back and user selection is based on 
either one or a combination of these components |4l- |9l: 

1) Channel Direction Information (CDI), e.g., beam index(BI), quantized channel index (QCI),..., etc. 

2) Channel Quality Information (CQI), e.g., SNR, SINR, channel norm, etc. 

3) User identity (ID) 

Another way to differentiate Feedback schemes according to whether their feed back data is analog 
or digital ifTTTl - lfl4l . While the forward link is opportunistic in nature, almost all feedback schemes 
mentioned above are non-opportunistic, and so, each user has a dedicated feedback channel. To see this, 
consider the random beamforming (RBF) scheme proposed by Sharif and Hassibi [4]. Here, the users are 
differentiated according to the channel direction (i.e. what beam direction is the channel mostly aligned 
with) and as a result, users feedback to Base Station (BS) the SINR corresponding to that direction only. 
So, each user feeds back one integer and one real number. In order to reduce the feedback load further 
Diaz et. al. in [7 ] propose a threshold based RBF. Here, instead of feeding back the SINR for the best 
beam for each receive antenna (one real plus one integer numbers), the user only transmits one bit to the 
BS, indicating whether or not the SINR on a pre-selected beam for any receive antenna is above a given 
threshold. The scheme is repeated for each beam. Since interaction or cooperation among the competing 
users is not allowed, hence defying opportunism, there is a linear increase in the feedback resources (or 
channels) with the number of users flU, 0. Even if thresholding is applied, there is no reduction in the 
number of feedback channels. This is because the channels are reserved even when users are not sending 
any feedback information. 

Recently, some works have started to consider opportunistic feedback schemes where feedback re- 
sources are shared and are opportunistically accessed by strong users, i.e. users whose CQI is above the 
given thresholds. Thus, in (8), Tang et. al. propose a feedback scheme with fixed number of feedback 
slots (channels) that are randomly accessed by strong users. In every slot, each strong user independently 
attempts to send back to the BS a data package containing its user identity (ID) with a probability. If two 
or more users feedback in the same slot, collision occurs and the feedback in that slot is discarded. When 
multiple users successfully feed back, the BS randomly selects one of the successful users. Although the 
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scheme requires only an integer feedback per slot, it is suboptimal as the user is selected randomly. The 
scheme was only proposed for the single-input single-output (SISO) case. In (9), Rajiv et. al. propose 
a similar feedback scheme based on random access slots for the MIMO case that requires only user 
identity feedback (i.e. an integer feedback per slot). Both of the two schemes above require accurate 
timing-synchronization to avoid collisions, which is difficult to achieve in practice. Moreover, feedback 
to the BS is successful if there are no collisions, i.e., only one user is attempting to feed back in a slot. In 
addition, the two schemes only work for digital feedback but not for analog. In all the feedback schemes 
discussed above, the feedback links were assumed ideal whereas the forward links were subject to both 
fading and noise. This asymmetry in the way the two links are treated is unrealisticO- 

In this paper, we consider a broadcast scenario where the forward and the feedback links are symmetric 
in that they are both i) non-ideal and ii) opportunistic or shared. Thus, both links undergo Rayleigh 
fading and are subject to additive Gaussian noise. Moreover, the channels in both links are shared and 
are opportunistic in the sense that feedback channels are dominated by strong users. Finally, the feedback 
links can be used for both analog and digital feedback. The paper proposes a generalized feedback model 
and compressive sensing (CS) [20] - E51 based opportunistic feedback protocols for feedback resource 
reduction. 

A. Our Feedback Approach 

To setup the stage, we summarize in what follows our feedback technique. This also helps in laying 
out our contributions and the organization of the paper. Just as in all existing feedback techniques, a 
number of channel directions or beams is first determined. For each direction, the number of feedback 
channels is fixed and strong users feedback their CQI information on all channels. In the analog feedback 
case, each strong user feeds back its CQI value whereas in the digital feedback case, each strong user 
feeds back "1" if its CQI is above a particular threshold and remains silent otherwise. This creates an 
undetermined system of equations in a sparse vector of users. We use the emerging compressive sensing 
technique to identify users who have fed back and to estimate the feedback CQI. Users with higher value 
CQI have a stronger chance of being recovered. The results obtained via compressive sensing are refined 
using least-squares and the strongest user is selected. As the feedback links are noisy, the BS backs off 
on the noisy CQI based on the variance of the noise. We obtain the optimum back off on the noisy CQI 
that maximizes the throughput. The scheme is repeated for each channel direction. Although we have 

'A notable exception is [35] and |36| which some how account for this by considering the effect of feedback error. 
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used SINR feedback, the proposed schemes can work with any kind of CQI (e.g. SNR). It is important 
to note that our scheme is less sensitive to timing-synchronization errors, as the scheme will be affected 
only if out of synchronization user is selected for a particular channel direction (the probability of which 
is low). 

The remainder of the paper is organized as follows. In Section JIJ a generalized feedback model is 
introduced. In Section HIT] we discuss the proposed feedback strategy. In Section [IVJ we present the sum- 
rate throughput obtained by the proposed schemes in the RBF case. In Section [V] performance evaluation 
of the proposed feedback schemes is presented. Feedback channel training is discussed in Section |VTJ 
followed by numerical results and conclusions in Sections IVIII and IVIIII respectively. 



A. Downlink Transmission Model 

We consider a single cell multi-antenna broadcast channel with p antennas at the base station (transmit- 
ter) and n users (receivers) each having one antenna^. The channel is described by a propagation matrix 
which is constant during the coherence interval and is known completely at the receiver. Let u G C pxl 
be the transmit symbol vector then the signal X{ received by the i-th user is given by 



where hj € C lxp is the channel gain vector between the transmitter and the user, and Z{ is the additive 
noise. The entries of hj and z% are i.i.d. complex Gaussian with zero mean and unit variance, CM (0, 1). 
Moreover, u satisfies unity average transmit power constraint E{u*u} = 1 and pi is the SNR of the i-th 
user. We consider a homogeneous network in which all users have the same SNR, i.e. pi = p = P/p, 
where P is the total power available at the transmitter. We also assume that the number of mobiles is 
greater than or equal to the number of transmit antennas, i.e., n > p, and that the BS selects p out of n 
users to transmit to. 

B. Generic Multi-Channel Feedback Model 

We present here a general model for the multiuser feedback channel with r feedback channels (possibly 
shared) among n users, in which users report channel quality information (CQI) to the base station in 

2 The broadcast model Q} is MISO. However, the feedback strategy in the paper applies to the general MIMO case. For in 
both cases, the users are feeding back a few number of coefficients, something that can be achieved by our feedback strategy. 



II. System Model 




i = 1, ... ,n 



(1) 
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order to exploit multiuser diversity. Let v £ c nxl be transmit feedback vector and let be the signal 
received via the i-th feedback channel described by 
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or equivalently 



Av + w 



(2) 



where r < n and aij represents the gain of the i-th channel for the j-th user. The feedback propagation 
matrix A is independent from the downlink channel and is assumed constant during the coherence interval 
and perfectly known at the BSq The entries of A are assumed to be zero mean unit variance complex 
i.i.d. complex Gaussian. The entries of the additive noise w are also complex i.i.d. Gaussian with zero 
mean and variance a 2 . 

We summarize in Table H how our feedback model © applies to the feedback models (opportunistic or 
not) suggested in literature flU- & Thus, in the non-opportunistic feedback model, each user is allocated 
its own feedback channel and the uplink channel matrix A becomes diagonal and of size n (equal to 
the number of users). For the opportunistic models proposed in |H1 by Tang et. al., the feedback channel 
matrix A becomes diagonal of size r x r, where r is the number of feedback slots and is less than n. 
The vector v represents feedback data in each slot, and when a collision in a particular slot takes place, 
the corresponding entry of v is not valid. The same model holds for |9l except that in this scheme r 
is not fixed but varies randomly. Also, r may not necessarily be less than n. In most of the proposed 
schemes, the additive noise w is set to zero. 

In this paper, we focus on the case when A is fading, the links are noisy, and r < n. Unlike the 
approach taken in ITT31 , fT6l , and |32l , which also consider noisy feedback, we consider a contention- 
based feedback protocol, which assigns independent multi-access contention channels for CQI reporting 
in which the feedback process should itself be a filter that selects strongest users. We assume that each 
single-antenna user is going to feedback the same information over r frequency bands shared with the 
other users. Thus, represents the gain of the j-th user in the i-th band (frequency feedback channels). 



3 We discuss in Section ITVl how this is done through feedback channel training. 
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III. Proposed Feedback Strategy 

Before we discuss the proposed feedback strategy, we present an important compressive sensing result 
that we will refer in our work. 

A. Compressive Sensing using LASSO 

We discuss an important result from a recent paper by Candes and Yaniv [25] for recovering v in a 
noisy setting (cf. ©) using the following LASSO estimate 

v = are; min II v — Avll? + aoilvll/, (3) 
veR" 11 J 11/2 11 11(1 v ' 

Assumptions: 

1) A has unit normed columns and obeys the coherence property. 

2) The noise is i.i.d. real Gaussian with variance a 2 , and v is sparse with support S. 

3) s < con/[||A|| 2 logn] for some positive numerical constant Co (where ||A|| is the operator norm of 
A), minjgs |v,| > 8a\/2 log n, and a is set to a = 2 y/2 log n. 

Under these assumptions, the LASSO estimate v obeys supp(v) = supp(v) and sgn(v) = sgn(v) with 
probability at least 1 — 2n _1 f-^=== + £ J — O ( n - 2 } og 2 ) . Moreover, once the support S is known, least 
squares can be used to estimate or refine V5 (non-zero entries of v) as = (A^A^) -1 A^y, where 
As denotes the sub-matrix formed by the columns {Aj\j G S}, indexed by the support S |[26l . 

B. General Strategy of Using Compressive Sensing for Feedback 

Any feedback scheme has two components, a direction component and a magnitude component. The 
transmitter usually has certain pre-determined directions for which it seeks user feedback. Thus, the BS 
announces that it is seeking feedback for a particular direction. The users whose channels lie at or are 
close to this direction, feedback their CQI (SNR, SINR, channel strength etc.). Now a limited number 
of users will feedback on the set of shared feedback channels according to equation ©. 

This makes the vector v in © sparse with sparsity level determined by the number of users who 
feedback. CS can now be used to recover the sparsity pattern of v. Moreover, the larger the value of 
particular CQI, the higher the chances of its recovery. Note that we need at least one strong user (i.e. 
s > 1) for each beam or direction in order to achieve full multiplexing gain which implies that small 
values of s are sufficient. To reduce the number of users who feedback s, we pursue a thresholding 
strategy where the user will feedback if its CQI is greater than a threshold ( to be determined. 
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Now consider a particular beam (CDI) (all beams will behave in an identical manner as the users are 
i.i.d. and the beams are equi-powered). Noting that the users' CQI are i.i.d., we can choose ( to produce 
a sparsity level s. This happen by requiring that 



F(Q = arg max ( '" )u s (l — u) r 
ue(o,i) \s,' 



n 



(4) 



where u = F(() or u is the complementary cumulative distribution function (CCDF) of CQI (SINR) 
defined a<H F(() = P[SINR > (] = ^feffi , C > 0. By taking the logarithm of both sides of ©, it is 
easy to show that the level needed to produce a sparsity level s is given by £ = F -1 (f-). 



C. Applying Candes and Yaniv result to our problem 

To apply Candes and Yaniv result to our problem, note first that the vector v in our case is real. So 
we increase the number of measurements by decomposing the observation y into its real and imaginary 
parts 



»(y) 
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where R(A) & 3(A) represents real and imaginary part of A, or, 

y = A v + w 



(5) 



The above model gives us the 2r x n real measurement matrix, and 2r x 1 real noise vector. To apply 
Candes and Yaniv result we still have to normalize the columns, thus 



y = Av + w 

where A = -4^ and v = yfrv. The outcome of the rearrangement gives us the following 



(6) 



1) The columns of A are i.i.d. Gaussian and normalized with ||A|| = y / n/2r. This implies that 

s < cq2t I 'logn or 2r > (l/co)slogn or 2r > cslogn. 

2) w are i.i.d. N (0, a 2 /2). 

3) mines' |vj| that can be recovered using LASSO is given by min^s |vj| > &a\J~^^ L . 

4 It was shown in |4| that the asymptotic behavior of the maximum of n i.i.d. random variables when n is sufficiently large 
converges in distribution to a random variable whose distribution function is exponential. The class of distribution functions we 
encounter in this paper are of the same type that were encountered in |4|. 
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D. Feedback Protocol for the Analog Feedback Case 

In the analog feedback scenario, users above threshold feed back their analog CQI value. For any 
given a, n and r, the threshold should be set such that mioses |vj| > 8a <J which implies that 
F~ l (s/n) > 8a<iJ^j^, or ( > 8a \ft^P L . Thus, only those value of s that satisfies the above equation 
are allowed. The other way of setting the threshold is to fix s, n and r and vary a (i.e. vary the SNR in 
the feedback channel ©) such that the above equation is satisfied. Note that s and r are also related by 
It > cs log n. 

The CS strategy then allows the BS to recover all users who transmitted their CQI with very high 
probability. This off course will be true provided that the number of users who feedback is less than or 
equal to s. Here, we assume that the probability a user is strongest for more than one beam is negligible 
as the number of users are relatively much larger than the number of beams. It has been shown in 
|4ll that this is a valid assumption under these conditions. The the proposed compressive sensing based 
opportunistic feedback protocol is shown in Table ITT] 

E. Feedback Protocol for the Digital Feedback Case 

The digital feedback is similar to analog feedback except that user feeds back "1" if its CQI for a 
particular beam is above a particular threshold. Otherwise, the user remains silent. Here, for a given n 



and r, a must be selected such that mhijgs |vj| > 8ay-^^ is satisfied which implies that a < 

From this inequality, it may first appear that s has no role to play here. However, note that s and r are 

related by 2r > cs log n. To increase the feedback granularity, we let the users compare their CQI to a 

set of thresholds, not just one. Thus, suppose that we want to set k thresholds £i < (2 < , < Ck such 

that the number of users whose CQI lie between the two consecutive thresholds Qi = [QXi+i) i s equal 
to s. Note that the last interval is 00) as (k+i = 00. Following our discussion in subsection IIII-BI 
we can set the lowermost threshold as F(d)n = sk, or, £1 = F _1 (^). Continuing in the same 
way, we get (2 = F~ Y ^ fc ~ 1 ^ , , Cfc = 1 (^)- The feedback procedure is as detailed in Table 

mil 

IV. Throughput in the RBF Case 

In this section, we present the sum-rate throughput achieved by the proposed schemes. Although we 
focus on RBF, the proposed scheme can be applied to other beamforming methods (e.g. ZFBF). 
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A. Throughput in the Analog Feedback Case 

The sum-rate throughput achieved in the RBF case with dedicated ideal feedback links is given by 



TZ^E 



V log 2 (l + max SINR ljm ) 

L — ' Ki<n 



,m=l 



pE 



log 2 (l+ max SINRi m ) 

Ki<n 



(V) 



where "sa" is used instead of "=" since there is a small probability that user i may be the strongest user 
for more than one beam. Using extreme value theory, it is shown in |4] that (O can be written as 



K = p log 2 (l + plog(n) - p(p - 2) log log(ra)) 
For the proposed scheme, the throughput in the analog case lZ a is given by 



n a = P E 



log 2 (l + max SINR im ) | A,B 

Ki<n 



E(AnB) 



(8) 



(9) 



where A is the event that CS is successful, and B is the event that for 1 < i < n, v, > (" is not a null 
set. As the events A and B are independent, we can write (© as 



TL a =pE 



log 2 (l+ max SINR iim ) | A,B 

Ki<n 



»(A)E(B) 



(10) 



-). To calculate E(B), note that the complement set 



As mentioned earlier, P(.A) = 1 — - . — 

v ' n \ ^2n log n 

B c is the event that the CQIs of all users are below the threshold (. This happens with probability 



P(B L ) = P( max Vi < C) = [1 - F(()] n = 1 - - 

l<i<n \ n 



s\ n 



Thus, 



n a = P E 



log 2 (l + max STNR iim ) 

l<j<71 



1-1 



1 



1 



n 



(11) 



n \\/27r logn / 

Due to noise, the fed back SINR might be higher than the actual one. This is due to the feedback 
channel noise, i.e. the feedback noise could add/subtract the original SINR value by a small value. As 
such, we need to estimate A correctly and subtract the effect of noise on the SINR. As the SINRs fed 
back by the users are transmitted as is and the feedback links are noisy, there is a need to back off the 
noisy received SINRs based on the noise variance. Specifically, the least square (LS) refined estimate is 
given by 

(AjAsr^y (12) 
(A^A 5 )- 1 A^(Av + w) (13) 
v 5 + e ' (14) 
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Here, the noise represented by e' is Gaussian as linear operations preserve the gaussian noise distribution; ■ 
A representative scalar equation of (fl4l) takes the form 

SINR' = SINR + e! (15) 

where SINR and SINR' stand for the actual and noisy SINR's respectively and where e' represents 
Gaussian noise. Now, if we decide to back off the received SINR by an amount A, then the back-off 
efficiency (rj) i.e. the probability that this backed off SINR is less than or equal to the actual SINR is 
given by 7/ = P[SINR' - A < SINR] = P[e' < A] = 1 - Q Thus, the effective throughput (with 

back-off on noisy SINR) can be written as 

K e// = (l ~ Q ( ^ H ft„ (16) 

= p (l-lA^^Vl-(l-infl-Q^W05-A) (17) 



n \\/27rlogn / / \ V nJ ) \ \<J e ' 
where f3 = 1 + p log(n) — p(p — 2) log log(n). We can now choose A to maximize R e ff- Differentiating 
Tleff w.r.t. A and setting it equal to yields 



A \ / /3 — A \ / A 



,2 



Q \^) + \^n\-Hr g{p - A) ^ <18) 

The value of A (which is a function of cr e <) that satisfies the above equation maximizes the effective 
throughput. It remains of course to calculate a e > which we undertake in subsection IV-BI further ahead. 

B. Throughput in the Digital Feedback Case 

The sum-rate throughput achieved by p beams in the multiple thresholds (k in number) based digital 
feedback case for RBF is given by 



n d ^ P E 



log 2 (l + max d) I A, B 

Ki<k 



where max Q is the lower limit of the CQI of the highest active threshold interval. Note that A and B are 

l<i<fc 

independent events, where A is the event that CS is successful and B is the event that maxK^fc Q is not 
a null set which basically means maxi<j< n Vj > £1 is not a null set. Thus, P(A) = 1 — ^ (^-j=== + i\ 
and P(B) (l - (l - f )") so that 

\ / 2 / 1 *\\ r 1 

(19) 



n ) ) V n V\/2vr logn n 



log 2 (l + max d 

Ki<k 



5 While the elements of e' are not necessarily independent, we will not use this fact in our backoff calculations. After all, these 
calculations are based on a worst case scenario. This will result in a larger backoff and a less optimistic achievable throughput. 
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We are left to evaluate E [log 2 (l + maxi<j<fc Q)] which can be derived analytically as (also see ETll ) 



E 



log 2 (l + max Ci) 

Kj<« 



log 2 (l + Ci)IP(selected user is in Qi)V{Qi) (20) 



i=i 



The probability of the threshold interval Qi is given by P(Qj) = [F(d + i) — F(d)], where F(() is 
the cumulative distribution function (CDF) of CQI (SINR) defined as F(() = 1 - ^q^#,C > 
J4[. Note that the probability that the selected user is in the threshold interval Qi is given as where 
V\ = P(j users other than the selected user are in Qi) = [F(& + i) — F(£j)p, and 
Vi = P((rt — j — 1) users lies below the interval Qi) = [F(Q)]^ n ~^~ 1 \ Substituting the values of V\ 
and V2 in (l20l . we find after straight forward manipulations 



P(selected user is in Q { ) - [F ^ i+1 ^ 1 f U ' !) ' 



n[F(Q +1 ) - F(C 
so that 



E 



log 2 (l + max Ci) 

Ki<k 



k 

- log 2 (l + Ci)([F(Ci+iT - [i ? (Ci)] n )and 



n 
i=l 



1= + £ )) - E lo g 2 (l + C,)([^(Cm)] n - [^(C,)] n ). 
log n n I J n f—f 

7 7 i=i 



V. Performance Evaluation 
We consider the following metrics for the performance evaluation of the proposed feedback schemes. 

A. Feedback Resources Reduction 

There is a significant reduction in number of feedback channels required for carrying feedback in- 
formation. The proposed schemes requires only 0(log(n)) feedback channels as opposed to n feedback 
channels required in the dedicated feedback casdj. More precisely, we have shown in subsection IIII-CI 
that the number of multiple access feedback channels should be at least r = |(slog(n)). 

6 We would like to stress here that we are reducing the resources dedicated to the user's feedback. In the independent feedback 
scenario that is usually pursued in the literature, the number of feedback channels is always n regardless of how small the 
feedback load is. 
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B. Feedback Noise Reduction in the Analog Feedback Case 

The other important benefit of this scheme is the feedback noise reduction (which eventually results in 
better throughput) in the analog feedback case. This is because the feedback data of each user is carried 
over all shared channels. This allows the base station to make use of the extra measurements available 
to reduce the effect of additive noise. This comes in contrast to dedicated channel feedback case where 
each user's feedback is carried over one channel only. In what follows, we study the error covariance 
that results from the LS refinement of the CS output. This will also allow us to determine the estimation 
error variance required in Subsection IIV-AI for SNR backoff in the analog feedback. 

1 ) Shared Feedback Channels: The error covariance matrix that results from the post-CS LS refine- 
ment of (fT2l) is given by 

R e = E As [R- 1 + A^, 1 A g ] _1 (21) 

where R v = E[v s v*] = er^I, and R w = E[w s w*] = cr^J and where S refers to the support indicating 
the users who feedback (obtained using CS). Thus, the S x S covariance matrix reads 



R e = a 2 v E As 



(22) 



Here p = is the signal to noise ratio. It is difficult to calculate this expectation in closed form. To 
go around this, note that what really matters to us are the diagonal values of R e which correspond to 
the effective estimated noise variance of the fed back SNR (see (fl9l)). It is not difficult to see that the 
maximum diagonal value of Re can be upper bounded by 



max a e (i) < a 2 E As 

Ki<s 



1 



(23) 



.1 ~\~ P^rnin{A-sA- s 

To calculate this expectation, we need to have the pdf of the minimum eigenvalue of the matrix A* A s . 
Now A s is a 2r x s matrix i.i.d. Gaussian matrix making A* A s a real central Wishart matrix W(s, 2r). It 
is difficult to obtain the pdf of \ m in and evaluate the expectation above for general J^, so we concentrate 
on the cases s = 1 and s = 2, and to understand the general trend we also consider the case of s/2r = f3, 
where r — > oo. Considering these cases makes sense given the fact that only a very low number of users 
are expected to feed back. We can show the following 

1) For s = 1, we have (Figure [7]) 



E 



r(l-r,l/2o) 



'Obtaining the pdf of \ m in for even the simple case of s = 2 proves to be very difficult to evaluate. 



December 30, 2010 



DRAFT 



13 



2) For s = 2, we have (Figure |7]) 



E 



l+pK 



1 



(2r-l)!(2r-2)! 



/o°° A S i r(2r + l)dA - 2 /» ^^r(2r, A)dA + / °° ^r(2r - 1, X)d\ 



poo A 2r - 1 e-* 



3) For s/2r = f5, with m — > oo, we have,(Figure [8) 



E 



i 



l+pA m! „ 



l+2pr(l- v / ?) 2 

2) Dedicated Feedback Channels: For the dedicated feedback case, the error covariance matrix is 
given by, 



ECM 



1 T 1 T 
— I + — 1 



(24) 



Thus, from (1241 . we conclude that the back off on the SINR is 0(^=) in the shared feedback channels 
case as opposed to 0(a w ) in the dedicated feedback channel case. 



C. Feedback Load Reduction 

In addition to the feedback resources reduction, there is a reduction in the amount of feedback. In 
RBF scheme with dedicated feedback channel Q, n real values and n integer values (nlog 2 p bits) are 
fed back, as there are n users in the system. 

1) Analog Feedback Case: The proposed CS based analog feedback scheme requires only pr real 
values to be fedback. This is because there are r shared channels and the scheme is repeated for each 
beam. Note that the feedback load reduction is more dominant in systems with large number of users, 
as r ~ 0(log(n)) and p is small. 

2) Digital Feedback Case: The proposed CS based digital feedback scheme requires only pkr bits 
to be fedback. This is because there are r shared channels and the scheme is repeated for each beam 
& threshold. Note that the feedback load reduction is more dominant in systems with large number of 
users, as r ~ 0(log(n)) and p & k are small. 

D. Trade-off in the Digital Feedback Case 

Given a budget of bits that can be fed back, using intuition, it was shown in ||3"2l that trade-off exists 
between the multi-user diversity and feedback accuracy. In our context, multi-user diversity is related to 
the number of shared channels r whereas feedback accuracy is related to the number of thresholds k, 
and so a similar trade-off may exist. The number of shared channels and thresholds must be chosen such 
that the throughput is maximized. This is explored using simulation in section IVIII 
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VI. Feedback Channel Training 

In the previous sections, we assumed that the channel A estimation is given to the system with the aid 
of a "genie" at no cost. In this section, we present how the feedback channel training can be accomplished 
and explore ways to reduce it. Here, we assume that (0 represents frequency feedback channels i.e., the 
entries of A, a^- represents the gain of the j-th user in the i-th frequency band. 

A. Channel Matrix is Full 

The optimal number of symbols required for channel training is at least equal to the number of transmit 
antennas [33]. So we need at least p training symbols for the downlink channel and at least n training 
symbol for the uplink channel (as there are n users each having one transmit antenna). In this section, 
we will discuss the uplink training for determining the entries of A and give an insight on the effect of 
error in A on compressive sensing. Note that training for each user in the uplink is performed one by 
one. So without loss of generality, we will discuss the training for user 1. 

Let us assume that training for user 1 is carried out using r symbols for the i-th frequency channel. 
Same procedure is repeated for all frequency channel. Mathematically we can write it as 

V% = y/f>ans T + Wi, i = l,...,r (25) 

where s T G C TXl is the training data (known at the BS), yi is the received signal at the BS, Wi is the 
gaussian noise and an is the feedback channel gain to be determined. 

The channel gains can be computed using the linear minimum-mean-square-error approach given in 
11331 . Thus, 

ail = Q + s T *s T ^ s T *yi, i = l,...,r (26) 

Thus, we can write 

an = an + ali, i = l,...,r (27) 

where a~n is the zero-mean channel estimation error. 

Note that training for all frequency channels can be done simultaneously. Thus, continuing in the same 
way for all users, we need nr symbols for the training. Also, it is important to note that we need one 
symbol for sending feedback data, so almost all of the uplink coherence time can be used for feedback 
training. Coherence time is typically of the order of few thousand symbols, so training would not be an 
issue for systems with moderate number of users. 



December 30, 2010 



DRAFT 



15 



In P4l . Tune et. al. showed that when compressive sensing is done using perfect A (i.e. with no 
estimation error), the number of channels satisfy the following 

Tperfect ~ ^tgil + f&) (28) 



where the denominator is the capacity of the channel with perfect channel gain a of a user. 

pa 2 
pa 2 +l 



With estimation error, denominator in the above equation modifies to log(l + ,, ). Note that this is 



obtained assuming a to be gaussian which is the worst case. Thus, 

r > lfl S("- fc + 1) . (29) 
" log(l + ^t) 

From the above equation, we can infer that estimation error results in the increase in the number of 
channels (for a fixed p) by a factor 

-!— = log(1 + pQ 2 2) . (30) 

T perfect log(l + ^j) 

The other way to mitigate the effect of estimation error is by increasing the SNR. Let us assume that 
Pperfect and p be the SNR in (l28l) and (|29l respectively. Then for a fixed number of channels, following 
equation must be satisfied 

1 ( P^ 2 \ nu 

Pperfect ~ a 2 \pa 2 + I J { ' 

A method for reducing the amount of feedback channel training time is discussed in the next subsection. 

B. Channel Matrix is Block Diagonal 

In order to reduce the feedback training time, we divide the users into groups with each group being 
allowed to feedback only on a set of feedback channels, thereby reducing the full channel matrix to a 
block diagonal one 



A 



BD 



A 



Compressive sensing is applied in the same way as discussed in Section Hill the only difference being 
that it is now applied on each block. Strong users in each block (or group) are found and the user 
corresponding to the maximum SINR among the strong users from all blocks is selected. As the users 
are i.i.d., so we divide the feedback resource equally among the k groups. Thus, training can now be 
performed for each block simultaneously. This approach reduces the feedback training time considerably, 
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e.g. if we divide the total number of users into two groups, then the training will require re/2 symbol 
time as opposed to re symbol time required for the case when the channel matrix is full. 

The flip side of this approach is that compressive sensing is now applied on the group of users instead 
of all users as one block. Thus, for same sparsity level s (overall), with block diagonalization, the number 
of feedback channels required is given below 

fA BD = /A.! H 1" /A fc 

= kf Al 




Note that from the above equation it may first appear that the number of channels have reduced as the 
quantity inside the logarithm is reduced by a factor of k, however, it is the other way round. This is 
because now c' has increased as the problem dimension (n) is reduced by a factor of k |[23l . Thus, there 
is a trade-off between the reduction in the amount of feedback training and the number of feedback 
channels. Also, note that there is now an additional constraint requiring s/k to be an integer. 

C. Non-fading Channels 

When the channels are non-fading i.e., the channel gains are constant (or 1), then each strong users 
multiples its CQI with a unique binary chip sequence (consisting of ±1 each with probability 0.5) of 
length equal to the number of shared feedback channels r and send it over the multiple access shared 
channels. There are two ways of assigning chip sequences to the users: pre-programmed in users' device 
or sending it over the air. If it is send over the air, then the training time (used in the case of fading 
channels) can be used to send unique binary chip sequences to all n users. Thus, A is r x n Bernoulli 
matrix and so CS can be applied as Bernoulli matrices are shown to satisfy the RIP [20]. 

VII. Numerical Results 

In this section, we present numerical results for CS-based feedback schemes by applying it in RBF 
context. We use p = 4 base station antennas, and n = 100 users. We set the threshold according to the 
sparsity level s, and use the maximum correlation technique (unless mentioned otherwise) for compressive 
sensing as this is much more computationally efficient than LASSO. Each point in the figures represents 
the sum-rate throughput achieved for shared number of channels determined by c & s. We use SNR = 10 
dB for both downlink and feedback link (unless stated otherwise) for calculating the sum-rate throughput. 
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A. Analog Feedback Case 

In Fig. [2j we present the sum-rate throughput with shared channel feedback in the analog feedback 
case. We use optimum back off on noisy SINRs in the analog feedback case. From this figure, we note 
that for small values of s the throughput is low. This is because the threshold works well for systems with 
large number of users but for systems with moderate number of users, we may have more or less number 
of users above the threshold than desired. So, if we set s low, then the probability that a beam has no 
strong user is relatively higher (resulting in a multiplexing loss) to the case when s is large. However, 
large values of s requires more feedback channels. Also, we see that the number of shared channels 
required to achieve the maximum possible throughput obtained in a noisy dedicated feedback scenario 
is 11 (corresponds to c/2 = 0.4 and s = 6). Also, it worth mentioning that the proposed scheme comes 
close to achieving the throughput obtained in a noiseless dedicated feedback scenario (dedicated feedback 
with ideal feedback links) due to feedback noise reduction. Note that 90% of throughput in noiseless 
dedicated feedback case is achieved by 19 shared channels (corresponds to c/2 = 0.8 and s = 5). 

In Fig. [3l we present results on block digitalization method proposed for reducing the feedback training 
time (section IVI-BI) . Here, we divide 100 users into two groups of 50 users and compressive sensing is 
applied on each group. It is clear from the figure that this method requires few more feedback channels. 
Also in this figure, we present result based on LASSO which shows that LASSO method performs 
marginally better than maximum correlation method. 

In Fig. |H we present the sum-rate throughput achieved by two-stage RBF in the analog feedback case 
when the feedback channels are noiseless. In the second-stage of the two-stage RBF, additional feedback 
information (beam gain information (BGI) [6]) is requested from the selected users only and Iterative 
Beam Power Control (IBPC) algorithm proposed in [6] is used for the re-distribution of the total power 
among the active beams (beams for which there are strong users) in an optimized manner. From the figure, 
we see that there is hardly any gain for two stage RBF with dedicated channel feedback, however, it is 
evident that for compressive sensing based opportunistic feedback protocol, two-stage RBF is effective 
even for moderate to large number of users. This is because if no user is strong for some beams, the 
system still suffers from the multiplexing loss but the power of those beam are distributed among the 
active beams in an optimized way. Also, note that there is no back-off required here as the feedback links 
are noiseless. With two-stage RBF, the number of shared feedback channels required is 15 (corresponds 
to c/2 = 0.8 and s = 4). 
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B. Digital Feedback Case 

For all digital feedback cases, we chose s = 1 (the minimum possible value) and set multiple thresholds 
as discussed in section IIII-Ei This is because for the proposed scheme, s = 1 will allow us to set the 
highest possible uppermost threshold thereby ensuring a higher throughput. 

In Fig. [51 we present the sum-rate throughput achieved with shared channel feedback in the digital 
feedback case. Is is evident form the figure that the proposed scheme in a noisy scenario achieves the 
throughput obtained in a noiseless dedicated feedback scenario (dedicated feedback with ideal feedback 
links). Also, we see that the throughput increases with the increase in the number of shared channels & 
thresholds. Taking the pessimistic view, we need only 10 feedback channels (corresponds to c/2 = 2 and 
s = 1). However, it is important to note that beyond a certain number of shared channels or thresholds, 
the throughput either becomes stagnant or increases marginally. 

In Fig. [6l we consider fixed budgets of p x kr bits that can be fed back. From the figure, we note 
that such a trade-off exists and for a given fixed budget there is an optimum number of thresholds and 
shared feedback channels that maximizes the throughput. 

VIII. Conclusions 

In this paper, a generic feedback channel model and compressive sensing based opportunistic feedback 
schemes are proposed. The proposed generic feedback channel model is shown to encompass all existing 
feedback channel models proposed in the literature. We have shown that the proposed analog & digital 
opportunistic feedback schemes achieves the same sum-rate throughput as that achieved by dedicated 
feedback schemes, but with feedback channels growing only logarithmically with number of users. Also, 
we derived an expression for the sum-rate throughput in the digital feedback case with multiple thresholds. 

In the analog feedback case (noisy scenario), it has also been shown that due to feedback noise 
reduction, the proposed scheme comes close to achieving the throughput obtained in the case of noiseless 
dedicated feedback. In the digital feedback case, it has also been shown that beyond a certain number 
of shared channels or thresholds, the throughput either becomes stagnant or increases marginally. Also, 
given a budget on the amount of bits that can be fedback, we have shown that there exist a trade-off 
between the number of shared channels and thresholds and therefore they must be chosen such that the 
throughput is maximized. 

Although the results presented here only show the performance of the the proposed schemes in the 
RBF context, the schemes can easily work with other beamforming methods. 
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Number of Shared Feedback Channels (r) c/2 



Figure 1 . Trace of Error Covariance Matrix for shared feedback 
channel (SFC) and dedicated feedback channel (DFC), n = 100 
and uplink SNR = 10 dB for different values of s. 



Figure 2. Analog Shared Channel Feedback: Throughput versus 
c/2 for RBF, p = 4, n = 100 and SNR = 10 dB (both downlink 
& feedback link) for different values of s. 
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Figure 3. Analog Shared Channel Feedback: Throughput versus 
c/2 for RBF, p = 4, n = 100 and SNR = 10 dB (both 
downlink & feedback link) for different methods viz. LASSO, 
maximum correlation, and maximum correlation with block 
diagonalization. 



Figure 4. Analog Shared Channel Feedback: Throughput versus 
c/2 for Single stage and Two stage RBF, p = 4, n = 100 and 
SNR = 10 dB & cxo for downlink and feedback link respectively. 
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Figure 5. Digital Shared Channel Feedback: Throughput versus 
c/2 for RBF, p = 4, n = 100 and SNR = 10 dB (both downlink 
& feedback link) for different values of k. 



Figure 6. Digital Shared Channel Feedback: Throughput versus 
k for RBF, p = 4, n = 100 and SNR = 10 dB (both downlink & 
feedback link) for different budgets of bits that are to be fedback. 
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Table I 

Generic Feedback Channel Model 



Table II 

Proposed compressive sensing based opportunistic feedback protocol for analog feedback. 



1) 


Threshold Determination: BS decides on thresholding level ( based on 


the sparsity level that can be 




recovered. 




2) 


User Feedback: Repeat the following steps for each beam. 






• CQI Determination: Each user determines its best beam (corresponding 


to the highest CQI value). 




• CQI Feedback: Each user feeds back its CQI if it is higher than £ on 


all shared channels. Otherwise, 




the user remains silent. 






• Compressive Sensing: BS finds the strong users using CS. 






• Least-squares estimation/refining: BS estimates or refines results obtained via CS using least-squares. 




. Optimum CQI Back off: BS backs off on the noisy CQI (SINR) based 


on the noise variance such that 




the throughput is maximized. 




3) 


User Selection: Select users and schedule them to beams. 
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Table III 

Proposed compressive sensing based opportunistic feedback protocol for the digital feedback. 



1) Threshold Determination: BS decides on thresholding levels £1, £2, (k based on the sparsity level that 
can be recovered. For each threshold interval Ci+i), repeat the User Feedback step. 

2) User Feedback: Repeat the following steps for each beam. 

• CQI Determination: Each user determines its best beam (corresponding to the highest CQI value). 

• CQI Feedback: Each user feeds back its CQI if it lies in threshold interval [Q, Ci+i) on all shared 
channels. Otherwise, the user remains silent. 

• Compressive Sensing: BS finds the strong users using Compressive Sensing. 

3) User Selection: For each beam, BS randomly selects one of strong users of the highest active threshold 
interval, where active threshold interval here means that there is at least one user sending feedback data in 
the interval. Here, CQI is the lower limit of the highest active threshold interval. 
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